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PROCEDE D'ESTJMATION DE-DISTANCE-POUR UN MOBILE SOUMIS A 

DES CONTRAINTES DYNAMIQUES DE PARCOURS 

L'invention concerne la navigation de terrain pour un mobile 
5 soumis a des contraintes de parcours variant au cours du temps, tel qu'un 
aeronef limite en taux de montee, la limite pouvant etre negative, et evoluant 
au-dessus d'une zone de terrain presentant des reliefs ou des obstacles 
menacants proches ou superieurs a son altitude de vol. 

Divers systemes ont ete developpes pour prevenir I'equipage d'un 
10 aeronef d'un risque de collision avec le sol. Certains, tels que les systemes 
TAWS (acronyme de I'expression anglo-saxonne : 'Terrain Awareness and 
Warning System"), font une prevision de trajectoire a court terme pour 
I'aeronef a partir des informations de vol (position, cap, orientation et 
amplitude du vecteur Vitesse) fournies par les equipements du bord, la 
15 placent en situation par rapport a une carte de la region survolee extraite 
d'une base de donnees d'elevation du terrain accessible du bord et emettent 
des alarmes a destination de I'equipage de I'aeronef a chaque fois que la 
trajectoire previsible a court terme entre en collision avec le sol. Ces 
systemes TAWS agrementent leurs alarmes de recommandations 
20 rudimentaires du genre 'Terrain Ahead, Pull up". Certains d'entre eux 
donnent egalement des informations sur le niveau de risque de collision que 
font encourir les reliefs et les obstacles environnant I'aeronef sous la forme 
d'une carte presentant les reliefs ou les obstacles du terrain survole en 
strates de couleurs differentes. Cependant, cette carte de risques de collision 
25 avec I'environnement tient uniquement compte des altitudes du relief 
relativement a la position du mobile et ne prend pas en compte de I'existence 
ou non d'une trajectoire realiste permettant de rejoindre les zones affichees. 

Pour satisfaire ce besoin de connaitre les points du terrain survole 
restant accessibles apres une manoeuvre d'evitement d'un relief ou d'un 
30 obstacle au sol, la_carte de risque de_ collision avec I'e.nvironnemeat dfilL 
afficher uniquement les zones pour lesquelles il existe un chemin possible 
depuis la position courante du mobile. La realisation d'un tel affichage passe 
par ('association d'une metrique a une carte du relief tiree d'une base de 
donnees d'elevation du terrain. 
35 Un © methode connue pour associer une metrique a une carte du 

relief tiree d'une base de donnees d'elevation du terrain a maillage regulier 
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deTa surface terrestre ou d'une partte-de celle=circorrsiste-a considerer la 
carte presentant le relief a partir de valeurs d'altitude figurant, avec les 
coordonnees geographiques : latitude et longitude des points de mesure, 
dans les elements de la base de donnees d'elevation du terrain comme une 

5 image dont les pixels sont les valeurs d'altitude des points de la base de 
donnees d'elevation du terrain illustrees dans la carte avec, pour 
coordonnees abscisse et ordonnee au sein de I'image, les coordonnees 
geographiques latitude et longitude de ces points figurant dans les elements 
de la base de donnees d'elevation du terrain et a faire appel a une 

10 transformee de distance operant par propagation pour estimer des distances 

au sein de cette image. 

Les transformees de distance operant par propagation egalement 
connues sous I'appellation "transformees distance de chanfrein" ("chamfer 
distance transform" ou "chamfer euclidean distance transform" en anglo- 

15 saxon) deduisent la distance d'un pixel dit pixel but par rapport a un autre 
pixel dit pixel source, des distances precedemment estimees pour les pixels 
de son voisinage, par un balayage des pixels de I'image. Le balayage permet 
d'estimer la distance d'un nouveau pixel but par rapport au pixel source par 
recherche du trajet de longueur minimale allant du nouveau pixel but au pixel 

20 source en passant par un pixel intermediate de son voisinage dont la 
distance a deja et6 estimee, la distance du nouveau pixel but a un pixel 
intermediaire de son voisinage dont la distance a deja ete estimee etant 
donnee par application d'un masque de voisinage communement appele 
masque de chanfrein. 

25 Une transformee de distance de ce genre a ete proposee en 1986 

par Gunilla Borgefors pour estimer des distances entre objets dans une 
image numerique, dans un article intitule :" Distance Transformation in Digital 
Images." et paru dans la revue "Computer Vision, Graphics and Image 
Processing", Vol. 34 pp. 344-378. Un des interets de ces transformees de 

30 distance par propagation est "de~reduFe la^complexite des calculs d'une 
estimation de distance en autorisant I'emploi de nombres entiers. 

Pour selectionner le trajet de longueur minimale donnant 
I'estimation distance, une transformee de distance par propagation doit tester 
tous les trajets possibles. Cette obligation se traduit par une contrainte de 

35 regularity imposee a I'ordre de balayage des pixels d'une image. G. 
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Borgefors proposeTpbur satisfaire cette conirairite de regularite7de"baiayer 
les pixels d'une image deux fois consecutivement, dans deux ordres inverses 
Tun de I'autre, qui sont soit I'ordre lexicographique, I'image etant analysee de 
gauche a droite ligne par ligne et de haut en bas, et I'ordre lexicographique 
5 inverse, soit I'ordre lexicographique transpose, I'image ayant subi une 
rotation de 90°, et I'ordre lexicographique transpose inverse. Elle propose 
egalement d'adopter un masque de chanfrein de dimensions 3x3 avec deux 
valeurs (3, 4) de distances de voisinage ou de dimensions 5x5 avec trois 
valeurs (5, 7, 11) de distances de voisinage. 

10 Les transformees de distance operant par propagation sont deja 

employees dans le domaine de la navigation terrain pour robot. Dans ce 
cadre, il est connu d'utiliser la transformee de distance de G. Borgefors avec 
une contrainte statique consistant a attribuer, de maniere autoritaire, une 
distance infinie a un point en analyse lorsqu'il apparaTt qu'il fait partie de 

15 reliefs ou d'obstacles a contourner repertories dans une memoire de zones 
interdites de franchissement, cela de maniere a eliminer, de I'ensemble des 
trajets testes lors d'une estimation de distance, ceux passant par les reliefs 
ou obstacles que le robot doit contourner. Cependant, une transformee de 
distance operant par propagation utilisee avec une contrainte statique dans 
20 le cadre de la navigation terrain pour robot, ne convient pas a la navigation 
terrain pour aeronef pour lequel la menace presentee par un relief ou un 
obstacle au sol depend du profil vertical de sa trajectoire. 

La presente invention a pour but un proced6 d'estimation des 
25 distances des points d'une carte extraite d'une base de donnees d'elevation 
du terrain par rapport a un point de reference employant une transformee de 
distance operant par propagation, avec une contrainte dynamique evolutive 
au cours du temps convenant a la navigation terrain pour un aeronef ayant 
une trajectoire a profil vertical impose. 



Elle a pour objet un procede d'estimation des distances des points 
d'une carte extraite d'une base de donnees d'elevation du terrain, pour un 
mobile soumis a des contraintes dynamiques lui interdisant certaines zones 
de la carte dites zones interdites de passage dont la configuration varie en 
35 fonction du temps de parcours du mobile. La base de donnees d'elevation du 
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~ terrain renferftr<run ensemble de-points reperes par une altitude^-une latitude- - 
et une longitude maillant le terrain devolution du mobile. Le precede met en 
oeuvre une transformee de distance operant par propagation sur I'image 
- constitute par les elements de la base de donnees d'elevation du terrain 

5 correspondent a la carte, disposes en lignes et colonnes par ordres de 
valeurs de longitude et de latitude. Cette transformee de distance estime les 
distances des differents points de I'image par rapport a un point source place 
a proximite du mobile en appliquant, par balayage, un masque de chanfrein 
aux differents points de I'image. L'estimation de distance d'un point, par 

10 application du masque de chanfrein a ce point dit point but s'effectue en 
repertoriant les differents trajets allant du point but au point source et passant 
par des points du voisinage du point but qui sont recouverts par le masque 
de chanfrein et dont les distances au point source ont ete prealablement 
estimees au cours du meme balayage, en determinant les longueurs des 

15 differents trajets repertories par sommation de la distance affectee au point 
de passage du voisinage et de sa distance au point but extraite du masque 
de chanfrein, en recherchant le trajet le plus court parmi les trajets 
repertories et en adoptant sa longueur comme estimation de la distance du 
point but. Initialement, en debut de balayage, une valeur de distance 

20 superieure a la plus grande distance mesurable sur I'image est attribute a 
tous les points de I'image sauf pour le point source, origine des mesures de 
distance, auquel est affectee une valeur de distance nulle. Le precede est 
remarquable en ce que les longueurs des trajets repertories, lors de 
I'application du masque de chanfrein a un point but, en vue de la recherche 

25 du trajet le plus court, sont traduites en temps de parcours pour le mobile et 
en ce que les trajets repertories dont les temps de parcours pour le mobile 
sont tels que le point but appartiendrait a une zone interdite de passage au 
moment oil le mobile I'atteindrait, sont exclus de la recherche du trajet le plus 
court. 

Avantageusement, lorsque le mobile est un aeronef ayant un profil 
vertical de vol a respecter determinant revolution de son altitude instantanee, 
on associe, aux longueurs des trajets repertoriees, les valeurs previsibles 
des altitudes instantanees qu'aurait I'aeronef en atteignant le point but par 
35 ces trajets tout en respectant le profil vertical de vol impose, et on elimine de 



WO 2005/031262 



PCT/EP2004/052078 



la recherche du trajef le plus court, les trajets repertories associes~a~des 
valeurs previsibles d'altitude inferieures ou egales a celle du point but 
donnee par la base de donnees d'elevation du terrain et augmentee d'une 
marge de protection. 

5 

Avantageusement, lorsque le mobile est un aeronef ayant un profit 
vertical de vol impose, l'estimation de distance operee par propagation sur 
I'image constitute des elements de la base de donnees d'elevation du terrain 
correspondant a la carte est doublee d'une estimation de I'altitude previsible 
10 pour I'aeronef au droit des differents points de I'image en supposant qu'il 
suive le trajet le plus court selectionne pour l'estimation de distance et qu'il 
respecte le profil vertical de vol impose. 

Avantageusement, lorsque le mobile est un aeronef a profil 
15 vertical de vol impose et que l'estimation de distance est doublee par une 
estimation de I'altitude previsible de I'aeronef, les altitudes des differents 
points de la carte sont soustraites des estimations des altitudes previsibles 
pour I'aeronef en ces points pour obtenir des ecarts par rapport au sol. 

20 Avantageusement, lorsque le mobile est un aeronef a profil 

vertical de vol impose et que l'estimation de distance est doublee par une 
estimation de I'altitude previsible de I'aeronef, les altitudes des differents 
points de la carte sont soustraites des estimations des altitudes previsibles 
pour I'aeronef en ces points pour obtenir des ecarts par rapport au sol 

25 affiches sur la carte en strates de couleurs. 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain correspondant a la carte, en plusieurs passes 
30 successives selon des ordres differents? 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain appartenant a la carte, en plusieurs passes 
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successives selon des ordres differehtset de maniere repette jusqu'a ce que 
les estimations de distance obtenues se stabilised. 

Avantageusement, la transformee de distance par propagation 
5 balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain correspondent a la carte, en plusieurs passes 
successives selon des ordres differents dont I'ordre lexicographique, I'ordre 
lexicographique inverse, I'ordre lexicographique transpose et I'ordre 
lexicographique transpose inverse. 

10 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation terrain correspondant a la carte, en une serie de quatre passes 
repetee jusqu'a stabilisation des estimations de distances : 
15 . une premiere passe effectuee ligne par ligne de haut en bas de 

I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
20 a droite de I'image, chaque colonne etant parcourue de haut 

en bas, et 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut. 

25 

Avantageusement, la transformee de distance par propagation 
balaie les pixels de I'image constitute des elements de la base de donnees 
d'elevation du terrain correspondant a la carte, en une serie de huit passes 
repetee jusqu'a stabilisation des estimations de distances. : 
30 . une premiere passe^ffectueeligne^par ligne de haut en bas de 

I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 
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~ unenroisieme passe effectueecolonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de haut 
en bas, 

- une quatrieme passe effectuee colonne par colonne de droite a 
5 gauche de I'image, chaque colonne etant parcourue de bas en 

haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une sixieme passe effectuee ligne par ligne de bas en haut de 
10 I'image, chaque ligne etant parcourue de gauche a droite, 

- une septieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 
15 Q droite d e I'image, chaque colonne etant parcourue de bas en 

haut. 

D'autres caracteristiques et avantages de I'invention ressortiront 
de la description ci-apres d'un mode de realisation donne a titre d'exemple. 
20 Cette description sera faite en regard du dessin dans lequel : 

- une figure 1 represente un exemple de masque de chanfrein, 

- des figures 2a et 2b montrent les cellules du masque de 
chanfrein illustre a la figure 1, qui sont utilisees dans une 
passe de balayage selon I'ordre lexicographique et dans une 

25 Passe de balayage selon I'ordre lexicographique inverse, 

- une figure 3 est un diagramme illustrant les principals etapes 
d'un precede, conforme a I'invention, pour estimer la distance 
d'un point en tenant compte d'une contrainte dynamique au 
cours de I'application d'un masque de chanfrein, 

& : une figure 4 est un diagramme Illustrant une variante du 

precede d'estimation de la distance d'un point montre a la 
figure 3, et 

- une figure 5 est un diagramme des principals etapes d'un 
precede, conforme a I'invention, d'estimation, par propagation, 

« des distances de I'ensemble des points d'une carte tenant 
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compte d'une conTrainte dyriamique et mettant en ceovrerun 
procede d'estimation de la distance d'un point tel que ceux 
montres aux figures 3 et 4. 

5 La distance d'entre deux points d'une surface est la longueur 

minimale de tous les parcours possibles sur la surface partant de Tun des 
points et aboutissant a I'autre. Dans une image formee de pixels repartis 
selon un maillage regulier de lignes, colonnes et diagonales, une 
transformee de distance par propagation estime la distance d'un pixel dit 

10 pixel "but" par rapport a un pixel dit pixel "source" en construisant 
progressivement, en partant du pixel source, le plus court trajet possible 
suivant le maillage des pixels et aboutissant au pixel but et en s'aidant des 
distances trouvees pour les pixels de I'image deja analyses et d'un tableau 
dit masque de chanfrein repertoriant les valeurs des distances entre un pixel 

15 et ses proches voisins. 

Comme montre a la figure 1, un masque de chanfrein se presente 
sous la forme d'un tableau avec une disposition de cases reproduisant le 
motif d'un pixel entoure de ses proches voisins. Au centre du motif, une case 
affectee de la valeur 0 repere le pixel pris pour origine des distances 

20 repertories dans le tableau. Autour de cette case centrale, s'agglomerent 
des cases peripheriques remplies de valeurs de distance non nulles et 
reprenant la disposition des pixels du voisinage d'un pixel suppose occuper 
la case centrale. La valeur de distance figurant dans une case peripherique 
est celle de la distance separant un pixel occupant la position de la case 

25 peripherique concernee, d'un pixel occupant la position de la case centrale. 
On remarque que les valeurs de distance se repartissent en cercles 
concentriques. Un premier cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places soit sur la ligne 
sqit sur la colonne du pixel de la case centrale sont affectees d'une valeur de 

30 distance D1 . Un deuxieme cercle de quatre cases correspondant aux quatre 
pixels les plus proches du pixel de la case centrale places en dehors de la 
ligne et de la colonne du pixel de la case centrale sont affectees d'une valeur 
de distance D2. Un troisieme cercle de huit cases correspondant aux huit 
pixels les plus proches du pixel de la case centrale places en dehors de la 
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ligne, de la colonne et des diagonals du pixel de'la case centraie sont 
affectees d'une valeur D3. 

Le masque de chanfrein peut couvrir un voisinage plus ou moins 
etendu du pixel de la case centraie en repertoriant les valeurs des distances 
d'un nombre plus ou moins important de cercles concentriques de pixels du 
voisinage. II peut etre reduit aux deux premiers cercles formes par les pixels 
du voisinage d'un pixel occupant la case centraie ou etre etendu au-dela des 
trois premiers cercles formes par les pixels du voisinage du pixel de la case 
centraie mais il est habituel de s'arreter a trois premiers cercles comme celui 
represents a la figure 1. Les valeurs des distances D1, D2, D3 qui 
correspondent a des distances euclidiennes sont exprimees dans une 
echelle autorisant I'emploi de nombres entiers au prix d'une certaine 
approximation. C'est ainsi que G. Borgefors donne a la distance d1 
correspondant a un echelon en absclsse x ou en ordonnee y la valeur 5, a la 
15 distance d2 correspondant a la rac ine de la somme des canes des echelons 
en abscisse et ordonnee Jx 2 +y 2 la valeur 7 qui est une approximation de 
5V2 , et a la distance d3 la valeur 1 1 qui est une approximation de 5^5 . 

La construction progressive du plus court trajet possible allant a 
un pixel but, en partant d'un pixel source et en suivant le maillage des pixels 
20 se fait par un balayage regulier des pixels de I'image au moyen du masque 
de chanfrein. Initialement, les pixels de I'image se voient affecter une valeur 
de distance infinie, en fait un nombre suffisamment eleve pour depasser 
toutes les valeurs des distances mesurables dans I'image, a I'exception du 
pixel source qui se voit affecter une valeur de distance nulle. Puis les valeurs 
25 initiates de distance affectees aux points but sont mises a jour au cours du 
balayage de I'image par le masque de chanfrein, une mise a jour consistent 
a remplacer une valeur de distance attribute a un point but, par une nouvelle 
valeur moindre resultant d'une estimation de distance faite a I'occasion d'une 

nouvelle applicatipn_d^ masque de chanfrein au point but considere. 

Une estimation de distance par application du masque de 
chanfrein a un pixel but consiste a repertorier tous les trajets allant de ce 
pixel but au pixel source et passant par un pixel du voisinage du pixel but 
dont la distance a deja ete estimee au cours du meme balayage, a 
rechercher parmi les trajets repertories, le ou les trajets les plus courts et a 
35 adopter la longueur du ou des trajets les plus courts comme estimation de 



30 



WO 2005/031262 



PCT/EP2004/052078 



10 

"distance. Cela se fait en placanTTe pixel Bunfcnron vearestimer la distance 
dans la case centrale du masque de chanfrein, en selectionnant les cases 
peripheriques du masque de chanfrein correspondant a des pixels du 
voisinage dont la distance vient d'etre mise a jour, en calculant les longueurs 

5 des trajets les plus courts reliant le pixel a mettre a jour au pixel source en 
passant par un des pixels selectionnes du voisinage, par addition de la 
valeur de distance affectee au pixel du voisinage concerne et de la valeur de 
distance donnee par le masque de chanfrein, et a adopter, comme 
estimation de distance, le minimum des valeurs de longueur de trajet 

10 obtenues et de I'ancienne valeur de distance affectee au pixel en cours 
d'analyse. 

L'ordre de balayage des pixels de I'image influe sur la fiabilite des 
estimations de distance et de leurs mises a jour car les trajets pris en compte 
en dependent. En fait, il est soumis a une contrainte de regularity qui fait que 

15 si les pixels de I'image sont reperes selon l'ordre lexicographique (pixels 
classes dans un ordre croissant ligne par ligne en partant du haut de I'image 
et en progressant vers le bas de I'image, et de gauche a droite au sein d'une 
ligne). et si un pixel p a ete analyse avant un pixel q alors un pixel p+x doit 
etre analyses avant le pixel q+x. Les ordres lexicographique (balayage des 

20 pixels de I'image ligne par ligne de haut en bas et, au sein d'une ligne, de 
gauche a droite), lexicographique inverse (balayage des pixels de I'image 
ligne par ligne de bas en haut et, au sein d'une ligne, de droite a gauche), 
lexicographique transpos6 (balayage des pixels de I'image colonne par 
colonne de gauche a droite et, au sein d'une colonne, de haut en bas), 

25 lexicographique transpose inverse (balayage des pixels par colonnes de 
droite a gauche et au sein d'une colonne de bas en haut) satisfont cette 
condition de regularity et plus generalement tous les balayages dans 
lesquels les lignes et colonnes, ou les diagonales sont balayees de droite a 
gauche ou de gauche a droite. G. Borgefors preconise un double balayage 

30^ des pixels de I'image, une fois dans I'ordTe lexicographique et une autre dans" 
l'ordre lexicographique inverse. 

La figure 2a montre, dans le cas d'une passe de balayage selon 
l'ordre lexicographique allant du coin superieur gauche au coin inferieur droit 
de I'image, les cases du masque de chanfrein de la figure 1 utilisees pour 

35 repertorier les trajets allant d'un pixel but place sur la case centrale (case 
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index6e-par 0)-acrpixel source en passanrpar un pixel du voisinage dont la 
distance a deja fait I'objet d'une estimation au cours du meme balayage. Ces 
cases sont au nombre de huit, disposees dans la partie superieure gauche 
du masque de chanfrein. II y a done huit trajets r6pertories pour la recherche 
5 du plus court dont la longueur est prise pour estimation de la distance. 

La figure 2b montre, dans le cas d'une passe de balayage selon 
I'ordre lexicographique inverse allant du coin inferieur droit au coin superieur 
gauche de I'image, les cases du masque de chanfrein de la figure 1 utilisees 
pour repertorier les trajets allant d'un pixel but place sur la case centrale 

10 (case indexee par 0) au pixel source en passant par un pixel du voisinage 
dont la distance a deja fait I'objet d'une estimation au cours du meme 
balayage. Ces cases sont complementaires de celles de la figure 2a. Elles 
sont egalement au nombre de huit mais disposees dans la partie inferieure 
droite du masque de chanfrein. II y a done encore huit trajets repertories pour 

15 la recherche du plus court dont la longueur est prise pour estimation de la 
distance. 

La transformee de distance par propagation dont le principe vient 
d'etre rappele sommairement a ete concue a Porigine pour Panalyse du 
positionnement d'objets dans une image mais elle n'a pas tarde a etre 

20 appliquee a I'estimation des distances sur une carte du relief extraite d'une 
base de donnees d'elevation du terrain a maillage regulier de la surface 
terrestre. En effet, une telle carte ne dispose pas explicitement d'une 
metrique puisqu'elle est tracee a partir des altitudes des points du maillage 
de la base de donnees d'elevation du terrain de la zone representee. Dans 

25 ce cadre, la transformee de distance par propagation est appliquee a une 
image dont les pixels sont les elements de la base de donnees d'elevation du 
terrain appartenant a la carte, e'est-a-dire, des valeurs d'altitude associees 
aux coordonnees geographiques latitude, longitude des noeuds du maillage 
ou elles ont ete mesurees, classes, comme sur la carte, par latitude et par 
~3o~ lonpude croissantes ou decroissantes seiorTur i tableau i deux dimensions 
de coordonnees latitude et longitude. 

Pour une navigation terrain de mobiles tels que des robots, la 
transformee de distance par propagation est utilisee pour estimer les 
distances des points de la carte d'un terrain devolution extraite d'une base 

35 de donnees d'elevation du terrain par rapport a la position du mobile ou une 
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position proche. Dans ce casTil est connu de tenircompte de -contraintes 

statiques constitutes par des zones de la carte infranchissables par ie 
mobile en raison de leurs configurations accidentees. Pour ce faire, un 
marqueur de zone interdite est associe aux elements de la base de donnees 

5 d'elevation du terrain figurant dans la carte. II signale, lorsqu'il est active, une 
zone infranchissable ou interdite et inhibe toute mise a jour autre qu'une 
initialisation, de I'estimation de distance faite par la transformee de distance 
par propagation pour I'element pixel considere. 

Dans le cas d'un aeronef, les zones infranchissables evoluent en 

10 fonction du profil vertical impose a sa trajectoire si bien qu'une estimation de 
distance sous contraintes statiques au moyen d'une transformee de distance 
par propagation n'est pas satisfaisante. 

On propose de prendre en compte, dans la definition des zones 
interdites de passage, I'altitude previsible de I'aeronef a chaque point but 

15 dont la distance est en cours d'estimation. Cette altitude previsible, qui 
depend bien evidemment du trajet emprunte, est celle de I'aeronef apres 
suivi du trajet adopte pour la mesure de distance. L'estimation de cette 
altitude previsible de I'aeronef en un point but, se fait par propagation au 
cours du balayage de I'image par le masque de chanfrein d'une maniere 

20 analogue a l'estimation de distance. Pour chaque trajet repertorie allant d'un 
point but au point source en passant par un point du voisinage du point but 
dont la distance au point source et I'altitude previsible de I'aeronef ont deja 
ete estimees au cours du meme balayage, I'altitude previsible de I'aeronef 
est deduite de la longueur du trajet et du profil vertical impose a la trajectoire 

25 de I'aeronef. Cette altitude previsible, estimee pour chaque trajet repertorie 
allant d'un point but dont la distance est en cours d'estimation a un point 
source place a proximite de la position de I'aeronef, est utilisee comme un 
critere de selection des trajets pris en compte dans I'estimation distance. Si 
elle est inferieure ou egale a I'altitude du point but figurant dans la base de 

30 do^nnees d'elevatiorTdu terrain majoree d'une marge (Je "securite, le Trajet 
repertorie auquel elle est associee est ecarte et ne participe pas a la 
selection du plus court trajet. Une fois la selection du plus court trajet 
effectuee, sa longueur est prise pour distance du point but et I'altitude 
previsible pour I'aeronef qui lui est associee est egalement retenue pour 

35 I'altitude de I'aeronef au point but. 
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La figure 3 illustreries-principales etapes du traitement effectue" 
lors de I'application du masque de chanfrein a un point but Py pour estimer 
sa distance pour un aeronef ayant un profil vertical de trajectoire impose. Le 
point but consider Py est place dans la case centrale du masque de 
5 chanfrein. Pour chaque point voisin P v qui entre dans les cases du masque 
de chanfrein et dont la distance a deja ete estimee au cours du meme 
balayage, le traitement consiste a : 

- lire la distance estimee D v du point voisin P v (etape 30), 

- lire I'altitude Ay du point but Py dans la base de donnees 
1 o d'elevation du terrain (etape 31 ), 

- lire le coefficient Cxy du masque de chanfrein correspondant a 
la case occupee par le point voisin P v (etape 32), 

- calculer la distance propagee D P correspondant a la somme de 
la distance estimee D v du point voisin P v et du coefficient Cxy 

15 affecte a la case du masque de chanfrein occupee par le point 

voisin P V: 

(etape 33), 

- calculer I'altitude previsible A P de I'aeronef apres 
20 franchissement de la distance D P directement a partir de la 

distance D P si le profil vertical impose a la trajectoire de 
I'aeronef est defini en fonction de la distance parcourue PV(D P ) 
et prend implicitement en compte le temps de parcours ou 
indirectement par I'intermediaire du temps de parcours si le 
25 P«*ofil vertical impose a la trajectoire de I'aeronef est defini par 

une vitesse de changement d'altitude (etape 34), 

- comparer I'altitude previsible Ap obtenue avec celle Ay du point 
but Pjj tiree de la base de donnees d'elevation du terrain 
augmentee d'une marge de securite A (etape 35), 

30 - elirniherla disTance~propagee D P si I'altitude prevFsTbleAVest ~ 

inferieure ou egale a celle Ay du point but Py tiree de la base 
de donnees d'elevation du terrain et majoree par la marge de 
securite A (etape 36), 

- si I'altitude previsible A P est superieure a celle Ay du point but 
35 Pij majoree par la marge de securite A, lire la distance Dy deja 
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affectee au point but consrdere-Pjj< etape 37) et la comparer a — 

la distance propagee Dp (etape 38), 

- eliminer la distance propagee D P si elle est superieure ou 
egale a la distance Dg deja affectee au point but considere Pjj, 

5 et 

- remplacer la distance D tii deja affectee au point but considere 
Pij, par la distance propagee Dp si cette derniere est inferieure 
(etape 39). 

La figure 4 illustre les principals etapes d'une variante du 

10 traitement effectue lors de I'application du masque de chanfrein a un point 
but Pjj pour estimer sa distance pour un aeronef ayant un profil vertical de 
trajectoire impose. 

Cette variante differe dans la maniere d'elaborer I'altitude 
previsible A P de I'aeronef. Elle suppose que I'altitude previsible pour I'aeronef 

15 en chaque point de la base de donnees d'elevation du terrain calculee en 
fonction du profil vertical impose a sa trajectoire et a partir de la longueur du 
trajet selectionne pour la mesure de distance n'est pas consideree comme 
une variable passagere, ce que permet le traitement decrit relativement a la 
figure 3, mais est memorisee, au meme titre que I'estimation de distance. 

20 Dans cette variante, les etapes 30, 31 de lectures de la distance estimee D v 
du point voisin P v et de I'altitude Ay du point but Pjj dans la base de donnees 
d'elevation du terrain sont completees par une etape 40 de lecture de 
I'altitude previsible Apv pour I'aeronef au point voisin P v , et le calcul de 
I'altitude previsible A P se fait (etape 34') par sommation de I'altitude previsible 

25 Apv au point voisin P v et de la variation d'altitude sur la distance separant le 
point voisin P v du point but Pjj due au profil vertical impose a la trajectoire de 
I'aeronef. 

La memorisation des altitudes previsibles pour I'aeronef lorsqu'il 
atteint les differents points de la carte qui lui sont accessibles permet 

Ito d'etablir, en leur soustrayalft le^s alfitu^es des^ points de la carte tirees de 1a 
base de donnees d'elevation du terrain, une carte des possibilites maximales 
de hauteurs de survol de I'aeronef representant les ecarts pr6visibles par 
rapport au terrain en strates de couleurs. Une telle carte facilite pour 
I'equipage de I'aeronef le choix d'une trajectoire realiste presentant la 

35 meilleure garde au sol. 
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Comme indique precedemment, ('estimation des distances des 
differents points de la carte se fait en appliquant un traitement par masque 
de chanfrein tel que ceux qui viennent d'etre decrits relativement aux figures 
3 et 4, a I'ensemble des pixels de I'image formee par les elements de la base 
5 de donnees d'elevation du terrain appartenant a la carte, pris 
successivement selon un balayage regulier comportant un minimum de deux 
passes realisees dans des ordres inverses. 

La figure 5 illustre les principales etapes d'un exemple de 
processus global permettant I'estimation des distances de I'ensemble des 
10 points d'une carte de relief pour un mobile soumis a des contraintes 
dynamiques. 

La premiere etape 50 du processus est une initialisation des 
distances affectees aux differents points de la carte consideres comme les 
pixels d'une image. Cette initialisation des distances consiste, comme 
15 indique precedemment, a attribuer une valeur de distance infinie, a tout le 
moins superieure a la plus grande distance mesurable sur la carte, pour tous 
les points de la carte consideres comme des points but, a I'exception d'un 
seul considere comme la source de toutes les distances auquel est attribue 
une valeur de distance nulle. Ce point source est choisi a proximite de la 
20 position instantanee du mobile sur la carte. 

Les etapes suivantes 51 a 54 sont des passes d'un balayage 
regulier au cours desquelles le masque de chanfrein est applique 
successivement et a plusieurs reprises a tous les points de la carte 
consideres comme les pixels d'une image, ('application du masque de 
25 chanfrein a un point de la carte donnant une estimation de la distance de ce 
point par rapport au point source, par execution d'un des traitements decrits 
relativement a la figure 3 ou a la figure 4. 

La premiere passe de balayage (etape 51) se fait dans I'ordre 
lexicographique, les pixels de I'imag e etant analyses ligne par ligne du haut 
30 vers le bas de I'image et de gauche a droite au sein d'une meme ligne. La 
deuxieme passe de balayage (etape 52) se fait dans I'ordre lexicographique 
inverse, les pixels de I'image etant toujours analyses ligne par ligne mais du 
bas vers le haut de I'image et de droite a gauche au sein d'une ligne. La 
troisieme passe de balayage (etape 53) se fait dans I'ordre lexicographique 
35 transpose, les pixels de I'image etant analyses colonne par colonne de la 
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gauche vers la droite deTimage et de haut en toasrau sein"Tl , une-meme- 
colonne. La quatrieme passe de balayage (etape 54) se fait dans I'ordre 
lexicographique transpose inverse, les pixels de I'image etant analyses 
colonne par colonne mais de la droite vers la gauche de I'image et de bas en 

5 haut au sein d'une meme colonne. 

Ces quatre passes ( etapes 51 a 54) sont repetees tant que 
I'image de distance obtenue change. Pour ce faire, le contenu de I'image de 
distance obtenu est memorise (etape 56) apres chaque serie de quatre 
passes ( etapes 51 a 54) et compare avec le contenu de I'image de distance 

10 obtenu a la serie precedente (etape 55), la boucle n'etant brisee que lorsque 
la comparaison montre que le contenu de I'image distance ne varie plus. 

En theorie, deux passes de balayage selon I'ordre lexicographique 
et I'ordre lexicographique inverse peuvent suffire* Cependant la presence de 
zones interdites de passage de forme concave peut provoquer, dans le 

15 phenomene de propagation des distances, des angles morts renfermant des 
pixels pour lequel I'application du masque de chanfrein ne donne pas 
d'estimation distance. Pour diminuer ce risque d'angle mort, il y a lieu de faire 
varier la direction du phenomene de propagation distance en faisant varier la 
direction du balayage d'ou le doublement des passes avec une transposition 

20 des ordres de balayage correspondant a une rotation de I'image de 90°. Pour 
une encore meilleure 6limination des angles morts, on peut proceder a des 
series de huit passes : 

- une premiere passe effectuee ligne par ligne de haut en bas de 
I'image, chaque ligne etant parcourue de gauche a droite, 

25 - une deuxieme passe effectuee ligne par ligne de bas en haut 

de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
a droite de I'image, chaque colonne etant parcourue de haut 
en bas, 

30 - une quatrieme passe effectuee coldrme^a?~c6lonhe deUfoite a 

gauche de I'image, chaque colonne etant parcourue de bas en 
haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de I'image, chaque ligne etant parcourue de droite a gauche, 
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- unesixieme passe effectuee ligne par ligne de bas~ eri haufde 
Pimage, chaque ligne etant parcourue de gauche a droite, 

- une septieme passe effectuee colonne par colonne de droite a 
gauche de l'image, chaque colonne etant parcourue de haut en 

5 bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 
a droite de l'image, chaque colonne etant parcourue de bas en 
haut. 

II est possible d'introduire dans la serie de passes de balayage 
10 d'autres types de passes de balayage deduits des passes precedentes en 
faisant jouer par les diagonales de l'image, les rdles joues precedemment 
par les lignes et colonnes de l'image. Cela revient a appliquer les passes de 
balayage decrites precedemment a une image tournee de 45°. De maniere 
generate, plus les passes d'une serie sont variees plus le risque d'angle mort 
15 diminue. 
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REVINDICATIONS 

1 . Precede d'estimation des distances des points d'une carte 

5 extraite d'une base de donnees d'elevation du terrain, pour un mobile soumis 
a des contraintes dynamiques lui interdisant certaines zones de la carte dites 
zones interdites de passage dont la configuration varie en fonction du temps 
de parcours du mobile ; la base de donnees d'elevation du terrain renfermant 
un ensemble de points reperes par une altitude, une latitude et une longitude 

10 maillant le terrain devolution du mobile ; ledit precede mettant en oeuvre une 
transformee de distance operant par propagation sur I'image constitute par 
les elements de la base de donnees d'elevation du terrain correspondant a la 
carte, disposes en lignes et colonnes par ordres de valeurs de longitude et 
de latitude ; la transformee de distance estimant les distances des differents 

15 points de I'image par rapport a un point source place a proximite du mobile, 
en appliquant, par balayage, un masque de chanfrein aux differents points de 
I'image ; I'estimation de distance d'un point, par application du masque de 
chanfrein a ce point dit point but s'effectuant en repertoriant les differents 
trajets allant du point but au point source et passant par des points du 

20 voisinage du point but qui sont couverts par le masque de chanfrein et dont 
les distances au point source ont ete prealablement estimees au cours du 
meme balayage, en determinant les longueurs des differents trajets 
repertories par sommation de la distance affectee au point de passage du 
voisinage et de sa distance au point but extraite du masque de chanfrein, en 

25 recherchant le trajet le plus court parmi les trajets repertories et en adoptant 
sa longueur comme estimation de la distance du point but ; une valeur de 
distance superieure a la plus grande distance mesurable sur I'image etant 
initialement attribute, en debut de balayage, a tous les points de I'image sauf 
au point source, origine des mesures de distance, auquel est affectee une 

30 valeur de distance nulle_ ;_ledit precede _etant caracterise en ce quejes 
longueurs des trajets repertories, lors de I'application du masque de 
chanfrein a un point but, en vue de la recherche du trajet le plus court, sont 
traduites en temps de parcours pour le mobile et en ce que les trajets 
repertories dont les temps de parcours pour le mobile sont tels que le point 

35 but appartiendrait a une zone interdite de passage au moment ou le mobile 
I'atteindrait, sont exclus de la recherche du trajet le plus court. 
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2. Precede selon la revendication 1, applique a im aeronef ayant 
un profil vertical de vol a respecter determinant revolution de son altitude 
instantanee, caracterise en ce que, aux longueurs des trajets repertoriees 

5 lors de ('application du masque de chanfrein a un point but, sont associees 
les valeurs previsibles des altitudes instantanees qu'aurait I'aeronef en 
atteignant le point but par ces trajets tout en respectant le profil vertical de 
vol impose, et en ce que les trajets repertories associes a des valeurs 
previsibles d'altitude inferieures ou egales a celle du point but donnee par la 
10 base de donnees d'elevation du terrain et augmentee d'une marge de 
protection sont exclus de la recherche du trajet le plus court,. 

3. Precede selon la revendication 2, caracterise en ce que 
I'estimation de distance operee par propagation sur I'image constituee des 

15 elements de la base de donnees d'elevation du terrain correspondant a la 
carte est doublee d'une estimation de I'altitude previsible pour I'aeronef au 
droit des differents points de I'image en supposant qu'il suive le trajet le plus 
court selectionne pour I'estimation de distance et qu'il respecte le profil 
vertical de vol impose. 

20 

4. Precede selon la revendication 3, caracterise en ce que les 
altitudes des differents points de la carte sont soustraites des estimations 
des altitudes previsibles pour I'aeronef en ces points pour obtenir des ecarts 
par rapport au sol. 

25 

5. Precede selon la revendication 4, caracterise en ce que les 
ecarts par rapport au sol sont affiches sur la carte en strates de couleurs. 

6. Precede selon la revendication 1, caracterise en ce que la 
30 transformee de distance par propagation" balaie les~ pixels de I'image 

constituee des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en plusieurs passes successives selon des ordres 
differents. 
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7. Precede selon la revendication 6, caracterise ^n ce que 1a 
transformee de distance par propagation balaie les pixels de I'image 
constituee des elements de la base de donnees d'elevation du terrain 
appartenant a la carte, en plusieurs passes successives selon des ordres 

5 differents et de maniere repetee jusqu'a ce que les estimations de distance 
obtenues se stabilisent. 

8. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 

10 constituee des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en plusieurs passes successives selon des ordres 
differents dont I'ordre lexicographique, I'ordre lexicographique inverse, I'ordre 
lexicographique transpose et I'ordre lexicographique transpose inverse. 

15 g. Procede selon la revendication 6, caracterise en ce que la 

transformee de distance par propagation balaie les pixels de I'image 
constituee des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en une serie de quatre passes repetee jusqu'a 
stabilisation des estimations de distances : 

20 - une premiere passe effectuee ligne par ligne de haut en bas de 

I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une troisieme passe effectuee colonne par colonne de gauche 
25 a droite de I'image, chaque colonne etant parcourue de haut 

en bas, et 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 
haut. 

30 

10. Procede selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constituee des elements de la base de donnees d'elevation du terrain 
correspondant a la carte, en une serie de huit passes repetee jusqu'a 
35 stabilisation des estimations de distances : 
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- une premiere-passe effectuee ligne par ligne de~h'aut en bas de~~ 
I'image, chaque ligne etant parcourue de gauche a droite, 

- une deuxieme passe effectuee ligne par ligne de bas en haut 
de I'image, chaque ligne etant parcourue de droite a gauche, 

5 - une troisieme passe effectuee colonne par colonne de gauche 

a droite de I'image, chaque colonne etant parcourue de haut 
en bas, 

- une quatrieme passe effectuee colonne par colonne de droite a 
gauche de I'image, chaque colonne etant parcourue de bas en 

10 haut, 

- une cinquieme passe effectuee ligne par ligne de haut en bas 
de I'image, chaque ligne etant parcourue de droite a gauche, 

- une sixieme passe effectuee ligne par ligne de bas en haut de 
I'image, chaque ligne etant parcourue de gauche a droite, 

15 - une septieme passe effectuee colonne par colonne de droite a 

gauche de Timage, chaque colonne etant parcourue de haut en 
bas, et 

- une huitieme passe effectuee colonne par colonne de gauche 
a droite de Timage, chaque colonne etant parcourue de bas en 

20 haut. 

11. ProcedS selon la revendication 6, caracterise en ce que la 
transformee de distance par propagation balaie les pixels de I'image 
constitute des elements de la base de donnees d'elevation du terrain 
25 appartenant a la carte, en plusieurs passes successives selon des ordres 
differents dont certains consistent en un balayage de I'image par diagonales, 
d'un bord a I'autre et, au sein d'une diagonale, d'une extremite a I'autre. 
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Calcul de la distance propage D P 
Dp = D v +0^ 
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